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Abstract. The Large Area Telescope (LAT) aboard the Fermi satellite allows us to study 
the high-energy 7-ray sky with unprecedented sensitivity. However, the origin of 31% of 
the detected 7-ray sources remains unknown. This population of unassociated 7-ray sources 
may contain new object classes, among them sources of photons from self-annihilating or 
decaying non-baryonic dark matter. Fermi-LAT might be capable to detect up to a few of 
these dark matter subhalos as faint and moderately extended 7-ray sources with a tempo- 
rally steady high-energy emission. After applying corresponding selection cuts to the second 
year Fermi catalog 2FGL, we investigate 13 candidate objects in more detail including their 
multi-wavelength properties in the radio, infrared, optical, UV, and X-ray bands. For the 7- 
ray band, we analyze both the 24-month and 42-month Fermi-LAT data sets. We probe the 
7-ray spectra for indications of a spectral cutoff, which singles out four sources of particular 
interest. We find all sources to be compatible with a point-source scenario. Multi- wavelength 
associations and, in particular, their infrared color-color data indicate no source to be com- 
patible with a dark matter origin, and we find the majority of the candidates to probably 
originate from faint, high-frequency peaked BL Lac type objects. We discuss possibilities to 
further investigate source candidates and future prospects to search for dark matter subhalos. 
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1 Introduction 

Current all-sky surveys of the high-energy (HE) 7-ray sky provide unprecedented sensitivity 
to disseminate the population of high-energy 7-ray emitters. Based upon 24 months of 
data recorded with the Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space 
Telescope satellite (Fermi) [1], the Fermi Collaboration recently published the second Fermi- 
LAT point-source catalog (2FGL) [2]. The 2FGL contains 1873 sources detected between 
100 MeV and 100 GeV (with a significance S > 4<r), whereof approximately one third (576 
sources) are lacking reliable association with sources detected in other wavelength bands. 
On the contrary, the majority (~1 000) of the 1 297 associated sources have been classified 
to most likely originate from active galactic nuclei (AGN), in particular BL Lacs and flat 
spectrum radio quasars [3]. 

While it seems plausible that most of the unassociated (high-latitude) 7-ray sources 
are expected to originate from faint AGN, this sample may also contain new classes of 7- 
ray emitting sources [4-9]. In particular, this includes sources potentially driven by self- 
annihilating (or decaying) dark matter (DM), i.e., DM subhalos (see [10] and references 
therein). In this context, it is also interesting to note that the intensity of the isotropic 
diffuse 7-ray background [11] cannot be fully accounted for by the properties of known 7- 
ray emitters (e.g., blazars) extrapolated below the confusion limit of Fermi-LAT [12, 13]. 
Unravelling the nature of Fermi unidentified sources therefore remains a crucial task to 
tackle searches and constrains on new HE phenomena, e.g., self-annihilating DM. 

Significant evidence for the existence of a so far undiscovered form of matter, so-called 
DM, have been provided by various different astrophysical observations via its gravitational 
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imprint [14-16]. Cold DM manifests itself on both cosmological and galactic scales, i.e., 
prevailing the formation of large scale structures down to accounting for galactic halos and 
their DM substructure (DM subhalos) [17, 18]. Observations indicate DM as an unknown 
non-baryonic type of heavy, electrical neutral, and color neutral particle, which is very weakly 
interacting with standard model particles. A promising scenario comprises DM to be con- 
stituted by weakly interacting massive particles (WIMPs) of Majorana type, with a mass 
between a few hundred GeV and several TeV. With an interaction strength at the order 
of weak interactions, thermal production of WIMPs in the early Universe can consistently 
explain the measured DM relic density [19]. Appropriate WIMP candidates naturally arise 
in beyond standard model theories, e.g., Supersymmetry [20]. 

The discovery of WIMPs is then possible via three complementary approaches, i.e., 
collider, direct, and indirect detection techniques: while WIMPs may be directly produced 
in colliders with sufficient center-of-mass energy, e.g., the Large Hadron Collider (LHC) 
[21, 22], underground low- noise experiments [23] are sensitive to their scattering signatures 
with heavy nuclei. From the astrophysical point of view, WIMPs may be indirectly detected 
through their self-annihilation (or decay) to standard model final states, eventually producing 
7 rays, charged light hadrons and leptons, and neutrinos [24]. 

In this paper, we present a search for DM subhalo candidates in the 2FGL catalog, 
following up on the 1FGL catalog search we conducted in [10] , henceforth called Paper I; for 
related studies, see e.g. [8, 25-29]. The paper is structured as follows: In section 2, the 7-ray 
properties of detectable DM subhalos are summarized. The catalog search for candidate 
7-ray sources and a study of their spectral and temporal properties are described in section 
3, including the investigation of multi- wavelength counterparts and an analysis of archival 
UV and X-ray data. The results of this study are summarized and discussed in section 4. 

2 Gamma-ray emission of DM subhalos 

In the hierarchical formation of structures, galactic DM halos are anticipated to host a 
large population of smaller subhalos (up to 10 16 ). Their mass spectrum dN/dM follows 
a power-law distribution over the mass range M between lO^-lO -3 M and ~ 1O 1O M : 
dN/dM ex M~ a , where a £ [1.9; 2.0], see [17, 18, 30]. While some of the massive subhalos 
are expected to host the Milky Way's luminous dwarf spheroidal satellite galaxies (dSphs), 
baryonic content of low-mass subhalos and even concentrated massive subhalos can be lacking, 
see [31, 32] and references therein. In Paper I, we have shown that up to two massive DM 
subhalos between 10 5 and 10 8 M could be detectable with Fermi-LAT within the first two 
mission years, assuming common models on the self-annihilation of heavy WIMPs, their 
density distribution in DM subhalos, and the distribution of the subhalos in the Galaxy. 
These objects at distances of O(kpc) would appear in the 7-ray sky as moderately extended 
(6*68 ~ 0.5°) 7-ray sources above 10 GeV. 1 The faint, temporally constant 7-ray flux at 
energies E between 10 and 100 GeV is anticipated at the detection level of Fermi-LAT, (/> p (10— 
100 GeV) ~ 10~ 10 cm -2 s , owing to the small self-annihilation cross section of WIMPs, 
(o~v) ~ 3x 10~ 26 cm 3 s _1 . The differential 7-ray flux approximately follows a hard power-law 
(index T < 1.5) with a distinct cutoff to the WIMP mass m x , where we assumed WIMPs 
of m x = 500 (150) GeV annihilating to heavy quarks, gauge bosons (e.g., bb, W + W~), or 

1 Within the angle #68, a fraction of 68% of the total 7-ray luminosity is emitted. In comparison with the 
point spread function of Fermi-LAT, #68 = 0.5° corresponds to about 4<tpsf, where ctpsf ~ 0.13° at 10 GeV 
(see http:/ /www. slac.stanford.edu/exp/glast/groups/canda/lat .Performance. htm). 
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to the leptons r + r~, see figure 2 in Paper I. We emphasize that the expected 7-ray flux 
may be even higher when including the possible presence of sub-substructures (enhancement 
by a factor of 2 to 3) [33, 34] and (rather model-dependent) photon contributions from 
final state radiation and virtual internal Bremsstrahlung [35-38]. The rather low WIMP 
velocity in bound subhalos might also lead to resonances in the self-annihilation cross section 
(Sommerfeld enhancement) [39, 40]. Potential secondary emission from energetic charged 
leptons eventually produced by WIMP annihilation in Galactic photon and magnetic fields is 
expected to be rather faint and diffuse, e.g. [41-43], and we therefore anticipate DM subhalos 
to be 7-ray sources which are not associated to sources detected in other wavelength bands 
at lower energies. 

3 DM subhalo candidates in the 2FGL 
3.1 Source selection 

The selection of candidate objects in the class of unidentified sources listed in the 2FGL 2 is 
based upon the properties of 7-ray emitting DM subhalos discussed above: 

(i) The sample was reduced to sources at high galactic latitudes \b\ > 20°, to avoid con- 
fusion with conventional Galactic sources and to reduce the impact of diffuse Galactic 
7-ray emission. 

(ii) Sources were selected for a steady 7-ray flux, requiring the cataloged variability param- 
eter var < 41.64, which corresponds to a probability of P s > 1% for the source to be 
steady. 

(iii) To select sources potentially driven by massive WIMPs, a detection above 10 GeV was 
required. In addition, most of the high-energy pulsars located at high Galactic latitude 
are eliminated by this energy cut. Spectrally, 7-ray pulsars may resemble DM subhalos, 
given their stable emission of characteristically hard 7-ray spectra (r < 2) with typical 
cutoff energies E c between 1 and 10 GeV [44-47]. 

(iv) Spectrally hard sources were selected constraining the index of the cataloged power-law 
fit with T < 2.0. 

All but the last cut have been already used in Paper I. Applying cuts (i) to (iv) to the 
unassociated sources listed in the 2FGL, 14 unassociated 7-ray sources remain. With the 
exception of 2FGL J2339.6— 0532, all other sources have a HE flux close to the detection level 
of Fermi-L AT . This is consistent with the expectation for candidate sources as discussed in 
section 2. We therefore discarded the outstandingly bright object 2FGL J2339.6— 0532. Table 
1 lists the final sample of 13 candidates sources together with their positional and spectral 
properties. 

The source 2FGL J0031. 0+0724 has been extensively studied in Paper I, where also 
2FGL J0143.6-5844 has been listed. An updated discussion of 2FGL J0031. 0+0724 is pre- 
sented below. 2FGL J2257.9— 3646 has been claimed as a DM subhalo candidate in [29], 
while very high-energy (E > 100 GeV) follow-up observations of 2FGL J2347. 2+0707 have 
been conducted with MAGIC [27]. All remaining sources are new candidates. 

2 Version v06 
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2FGL name 


l,b 
fdegl 

L oj 


<7rs / (705 

[arcmin] 


5 

L J 


r 


P (1O-1OO GeV) 
[lO" 10 cm^s- 1 ] 


Ss 

L J 


J0031.0+0724 lst 


114.095,-55.108 


4.4/7.2 


4.4 


1.9 + 0.3 


0.7 + 0.3 


4.5 


J0116. 6-6153 


297.749,-54.986 


3.7/6.0 


5.5 


1.6 + 0.2 


0.7 + 0.3 


4.9 


J0143.6-5844 lst 


290.468,-57.102 


2.4/3.8 


14.2 


1.7 + 0.1 


2.2 + 0.6 


9.5 


J0305.0-1602 lst 


200.151,-57.146 


4.5/7.3 


5.3 


1.5 ± 0.2 


0.9 + 0.4 


4.4 


J0312. 8+2013 


162.507,-31.569 


3.7/6.0 


4.4 


1.7 + 0.2 


0.7 + 0.3 


4.5 


J0338. 2+1306 


173.471,-32.929 


3.9/6.3 


5.8 


1.5 ± 0.2 


1.1 + 0.5 


5.1 


J0438.0-7331 


286.088,-35.168 


4.1/6.6 


6.1 


1.4 + 0.2 


0.8 + 0.4 


5.0 


J0737.5-8246 


295.086,-25.467 


3.7/6.0 


4.4 


1.3 + 0.3 


1.0 + 0.4 


5.2 


J1223.3+7954 


124.470,+37.134 


3.6/5.8 


4.2 


1.4 + 0.3 


0.6 + 0.3 


4.6 


J1347.0-2956 


317.047,+31.398 


4.1/6.6 


5.0 


1.4 + 0.3 


1.0 + 0.5 


4.2 


J1410.4+7411 


115.839,+41.825 


2.9/4.7 


9.1 


1.9 + 0.1 


0.7 + 0.3 


4.7 


J2257.9-3646 lst 


3.899,-64.186 


5.0/8.2 


5.3 


1.9 + 0.2 


0.8 + 0.4 


4.1 


J2347.2+0707 lst 


96.214,-52.385 


3.7/5.9 


7.2 


2.0 + 0.2 


0.8 + 0.4 


4.1 



Table 1: DM subhalo candidates in the 2FGL catalog. The first column lists the 2FGL 
name, where the index "1st" flags sources which have already been listed in the 1FGL (i.e., 
1FGL J0030.7+0724, 1FGL J0143.9-5845, 1FGL J0305.2-1601, 1FGL J2257.9-3643, and 
1FGL J2347. 3+0710). For each source, the position is given in galactic coordinates (l,b), 
together with the positional uncertainty 068(95) [68% (95%) c.l., semi- major axis], detection 
significance S in Gaussian sigma, the power-law index T, and the photon flux in the 10- 
100 GeV band. The last column S5 lists the significance of the 10-100 GeV detection (in 
Gaussian sigma). 



3.2 Fermi-L AT data 

The data analysis of each object in table 1 was based on data recorded with the Fermi- 
LAT in the first 24 as well as 42 months 3 of the mission. 4 We chose the same analysis 
framework and recommended options that were used for the 2FGL (based upon 24 months of 
data), with the exception that the considered energy range was extended to cover 100 MeV to 
300 GeV. The data analysis was performed with the public version of the Fermi Science Tools 
(v9r23pl, release date 06 October 2011) using data of Pass- 7 event reconstruction along with 
the P7_V6 instrument response functions. 5 All events passing the SOURCE event class were 
considered. Events were filtered for a maximum zenith angle of 100° (to eliminate contamina- 
tion from the Earth's limb), a maximum rocking-angle of 52°, and the recommended quality 
filters DATA_QUAL == 1 and LAT_CONFIG == 1 were applied. The (quadratic) region- 
of-interest (Rol) was centered on the nominal 2FGL position of the source-of-interest (Sol) 
with a size of 20° x 20°. The overall spectral fit was performed using the binned likelihood 
method (with 10 energy-bins per decade; optimizer NEWMINUIT, requiring an absolute fit 
tolerance of 10~ 3 ), where the source model contained all 2FGL sources within the Rol, along 
with fixed cataloged positional and spectral parameters. Including the Sol, the normaliza- 
tions and indices V of the default model of a power-law spectrum [d(j)/dE = (fro (E / Eq)~ v ] 

3 The 42-month data set covers the time period between the beginning of August 2008 (239557419 MET) 
up to the beginning of February 2012 (350063020 MET). 

4 The LAT data is publicly available at http://fermi.gsfc.nasa.gov/ssc/data/. 
5 See http: //fermi. gsfc.nasa.gov / ssc/data/analysis / . 
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of the innermost six sources were kept free, while the energy scale Eq was fixed to the 
cataloged pivot energy. We used the latest publicly available models for the Galactic fore- 
ground {gaL2yearp7v6.vO.fits) and isotropic background emission (iso_p7v6source.txt) . The 
normalization and corrective power-law index of the Galactic foreground template and the 
normalization of the isotropic background template were left free. In detail, the analysis was 
performed with the tools gtselect, gtmktime, gtbin, gtltcube, gtexpcube2, gtsrcmaps, and gtlike. 
For each source, we checked that our overall fit reproduces the cataloged data sufficiently. 

The analysis chain for the 42-month data set allowed for a possible change of the posi- 
tional coordinates of each Sol. We used gtfindsrc to refit the position. The refined uncertainty 
contour (at 95% confidence level) was computed from the two-dimensional likelihood function 
£(RA, Dec), requiring 2A(log£) = 6.18 (2 degrees of freedom). 

3.3 Spectral analysis 

The energy spectra of the candidate sources have to be consistent with a spectrum gener- 
ated by self-annihilating WIMPs. For each candidate source, we carried out a statistical 
hypotheses test based upon a likelihood ratio. The null hypothesis (H p \) of the Sol to follow 
a conventional power-law spectrum was tested against the hypothesis (H exp ) of a spectrum 
generated by WIMP annihilation. The likelihood ratio defines a test statistic 



where C{H) denotes the total likelihood for the corresponding Rol, fitted assuming the Sol 
to follow the spectral hypothesis H (cf. section 3. 2). 6 

As benchmark models, we probed for WIMP annihilation to heavy quarks (e.g., bb), 
gauge bosons (e.g., VF + VF~), and to the leptons r + r~, which lead to a considerably harder 
7-ray spectrum. For our purpose, the differential 7-ray spectra resulting from annihilation 
of supersymmetric neutralinos to the afore mentioned final states, see [48, 49], can be ap- 
proximated with a power-law spectrum modified by an exponential cutoff [50], dN^/dx = 
Nq x~ T exp (— px), where x = E/m x . This simple parametrization originally introduced for 
gauge boson final states [50] also provides a reasonable fit to the r + r~ final states [49] with 
different values of Nq, T, and p. The used parameters are listed in table 2. This approach sim- 
plifies considerably the treatment and interpretation of the fit procedure. Given the small 
number statistics of the faint sources, more subtle differences in the final state spectrum 
cannot be resolved. 

Therefore, the hypothesis -ffexp was a power-law spectrum with exponential cutoff, 
d(j)/dE = 4>q (E/Eo)~ r exp (—E/E c ). The WIMP mass is then connected to the cutoff energy 
via m x = pE c , and the normalization is <pQ = AnTnE -1 Eq^ . The index T was fixed to the 
values given in table 2, while the energy scale was set to Eq = 1 GeV. As discussed in section 
3.1, we note that this hypothesis also probes for 7-ray pulsars which might contaminate our 
sample, even though we expect most of them to be eliminated by the energy cut. 

Since T was kept fixed for H exp , the null and alternative hypotheses are not nested. This 
implies that the test statistic TS ex p does not necessarily follow the theorems of Wilks [51] or 
Chernoff [52], i.e., TS ex p is not drawn from a chi-square distribution in the null hypothesis. 

6 Note that the spectral models (usually power laws) of all other sources in the Rol are kept fixed. Therefore, 
another possibility would be to use the test statistic TS assigned by the spectral likelihood fit (gtlike) to 
measure the source's detection significance, S oc VTS. We checked that both methods give similar results, as 
expected. 




(3.1) 
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Channel 


N 


r 


P 


Remarks 


heavy quarks, 


0.73 


1.5 


7.8 


[50] 


gauge bosons 








r + r - 


5.28 


0.35 


4.6 





Table 2: Fit parameters to approximate the differential photon spectra originating from 
final-state fragmentation of self-annihilating neutralinos by power laws with exponential cut- 
offs. For annihilation to heavy quarks and gauge bosons, we use the spectral parametrization 
from [50], while the photon yield from annihilation in t + t~ is approximated from [49]. 



It also implies that TS exp can have both negative as well as positive values: the hypothesis 
-ffexp is disfavored if TS exp <C and favored if TS exp S> 0. Since the distribution of the test 
statistic is a priori not known, the significance of this test had to be calculated with Monte 
Carlo simulations. The methods are described in appendix A. Based upon the simulations, we 
find that for the index T = 1.5 (0.35), a significance of 2a corresponds to TS cxp = —6 (—20) 
and TS exp = 2(2), respectively, while the 3a contour is given by TS cxp = —25 (—35) and 
TS CX p = 4(7); table 8 lists further significance and their corresponding TS exp values. 

The best-fit parameters for a fixed T = 1.5 and T = 0.35 are summarized in table 3 
for the 24-month data and in table 4 for the 42-month data. Although the fitted spectrum 
depends on just two free parameters, the statistical errors of the normalization and cutoff 
energy remain comparably large, owing to the sample of very faint sources studied. 7 For the 
24-month data, we find that a power-law with exponential cutoff spectrum is favored for the 
sources 2FGL J0305. 0-1602 and 2FGL J0338. 2+1306, with a significance of 2.4a and 3.3a, 
respectively. After 42 months, for both sources the significance of this initial indication has 
decreased. However, an exponential cutoff is now favored for the sources 2FGL J0143.6— 5844 
and 2FGL J1410. 4+7411, with significances of ~3<r. Vice versa, note that for no source we 
find such a spectrum to be disfavored by the data (i.e., by more than 3a). 

The initial indications for an exponential cutoff in the 24-month or 42-month data 
sets motivate a closer inspection of the sources 2FGL J0305. 0—1602, 2FGL J0338. 2+1306, 
2FGL J0143. 6-5844, and 2FGL J1410. 4+7411. A comprehensive discussion of each source 
candidate is given in section 3.6. 

3.4 Variability and angular extent 

A 7-ray signal from a DM subhalo is expected to be constant in time and may be resolved as 
an extended source. Therefore, the temporal and spatial distributions of high-energy photons 
were tested for compatibility with a constant flux and the hypothesis of angular extent. The 
method used is based upon our previous work in Paper I. Most of the source candidates 
have been detected exclusively in the upper energy bins with a comparably low background 
contamination. Therefore, the tests were applied to the high-energy photon distribution be- 
tween 3-300 GeV and the inclusive interval 10-300 GeV, respectively. As motivated in section 
2 (for details see Paper I), high-energy photons within a circular region of radius 0.5° around 
the nominal 2FGL position were examined. Due to the low background contamination, this 
sample is dominated by signal events (see below). 

7 Note that the errors quoted in tables 3 and 4 are correlated, given that the energy scale Eo was fixed to 
1 GeV. Choosing the decorrelation energy, in principle, reduces the statistical errors. 
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24 months 


r 

0o/lO- n 
[(cm 2 s GeV)- 1 ] 


= 1.5 

E c 
[GeV] 


rpo 

J- ^exp 


r = 

0o/lO- n 
[(cm 2 s GeV)- 1 ] 


0.35 

E c 
[GeV] 


J- ^cxp 


2b (jL name 


J0031.0+0724 


17.1 ±8.0 


67 ±84 


-2.2 


1.3 ± 1.1 


21 ± 12 


-5.8 


J0116. 6-6153 


20.0 ±9.0 


34 ±33 


0.9 


4.7 ±4.2 


9±5 


-0.4 


J0143.6-5844 


70.5 ± 13.2 


46 ±25 


1.4 


16.8 ±5.5 


10 ±2 


-17.2 


J0305.0-1602 


21.1 ±7.5 


70 ±78 


1.4 


4.8 ±2.7 


12 ±5 


4.1 


J0312.8+2013 


29.9 ± 11.7 


46 ±40 


-0.6 


2.8 ±2.1 


16 ±8 


-5.2 


J0338. 2+1306 


48.0 ± 13.5 


40 ±28 


4.2 


12.6 ±5.9 


9±3 


6.7 


J0438.0-7331 


24.5 ±9.4 


105 ± 134 


0.9 


3.8 ±2.5 


17 ±8 


1.1 


J0737.5-8246 


31. 6± 11.8 


66 ±68 


1.3 


5.0 ±2.8 


14 ±6 


2.6 


J1223.3+7954 


8.1 ±6.0 


186 ± 502 


0.1 


0.3 ±0.4 


41 ±35 


-0.2 


J1347.0-2956 


26.1 ±8.6 


204 ± 367 


0.2 


2.0 ± 1.2 


31 ±17 


-2.6 


J1410.4+7411 


40.4 ± 10.7 


30 ± 18 


1.2 


12.4 ±5.3 


7±2 


-3.7 


J2257.9-3646 


33.9 ± 12.6 


17 ± 12 


1.3 


9.4 ±7.0 


6±3 


-4.6 


J2347.2+0707 


36.6 ± 10.8 


61 ±49 


-2.5 


4.7 ±2.6 


15 ±6 


-12.1 



Table 3: Best-fit parameters for a power law with exponential cutoff spectrum, fitting the 
24-month data set between 0.1 and 300 GeV. The index was fixed to T = 1.5 or T = 0.35, 
respectively. The table lists the normalization 4>q and cutoff energy E c . The column TS cxp 
gives the likelihood ratio for the comparison with a pure power-law fit. See the text for more 
details on the interpretation of TS exp . 



42 months 


r 

^o/10- n 
[(cm 2 s GeV)" 1 ] 


= 1.5 

E c 
[GeV] 


TS CX p 


r = 

0o/lO- n 
[(cm 2 s GeV)- 1 ] 


0.35 

E c 
[GeV] 


TS CX p 


2FGL name 


J0031.0±0724 


17.5 ±7.7 


38 ±35 


-6.9 


1.1 ± 1.0 


18 ± 11 


-13.9 


J0116. 6-6153 


27.1 ±7.3 


44 ±30 


1.2 


5.8 ±3.0 


10 ±4 


-7.1 


J0143.6-5844 


53.6 ±8.7 


53 ±26 


3.9 


13.4 ±3.9 


10 ±2 


-17.8 


J0305.0-1602 




E c — > oo 




2.9 ± 1.7 


12 ±5 


2.3 


J0312.8±2013 


23.4 ±8.2 


42 ±32 


0.5 


3.0 ± 1.8 


13 ±5 


-3.0 


J0338.2±1306 


32.3 ±7.6 


153 ± 177 


1.1 


5.8 ±2.5 


14 ±5 


-2.7 


J0438.0-7331 


20.5 ±7.7 


54 ±48 


0.6 


2.8 ±1.9 


14 ±6 


-3.6 


J0737.5-8246 


29.1 ±9.4 


47 ±38 


2.2 


6.4 ±3.5 


10 ±4 


2.3 


J1223.3±7954 


7.8 ±4.9 


103 ± 179 


-0.3 


0.4 ±0.4 


31 ±23 


-2.1 


J1347.0-2956 


27.9 ±6.8 


212 ±285 


0.3 


2.1 ±0.9 


31 ± 13 


-7.5 


J1410.4+7411 


42.0 ±8.0 


27 ± 12 


4.3 


15.3 ±4.7 


7± 1 


-3.3 


J2257.9-3646 


19.5 ±7.9 


21 ± 15 


0.8 


3.3 ±2.4 


9±4 


-2.8 


J2347.2+0707 


46.9 ±9.5 


49 ±27 


-0.1 


12.0 ±4.7 


8±2 


-14.6 



Table 4: Best-fit parameters for a power law with exponential cutoff spectrum, fitting the 
42-month data set between 0.1 and 300 GeV. The index was fixed to T = 1.5 or T = 0.35, 
respectively. The table lists the normalization 0o and cutoff energy E c . The column TS exp 
gives the likelihood ratio for the comparison with a pure power-law fit. See the text for more 
details on the interpretation of TS exp . 



For the 24-month data, the upper part of table 5 lists the number of photons detected 
from each source between 3-10 GeV and 10-300 GeV, respectively, together with the energy 
and event class of the photon with highest energy. The lower part lists the same quantities 
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2b (jL name 


Number of photons (r < 0.5°) 


^(10-300 GeV) 


E 


Evcl 


6 — lu (je v 


±U — oUU (ore V 


gal/iso/E 


[GeV] 


J0031.0+0724 







0.3/0.6/0.9 


44 ±4 


4 


J0116.6-6153 




E 



0.2/0.7/0.9 


2S±\ 


4 


J0143.6-5844 


1 7 
1 1 


1 E. 


0.2/0.7/0.9 


45 ±3 


4 


J0305.0-1602 


1 1 
1 1 


D 


0.3/0.6/0.9 


39+? 


4 


J0312.8+2013 




E 




0.9/0.7/1.6 


35±| 


4 


J0338.2+1306 


27 


u 


1.3/0.6/1.9 


29 ±2 


4 


J0438.0-7331 


14 


6 


0.6/0.8/1.4 


56 ±3 


4 


J0737.5-8246 




7 


0.9/0.6/1.5 


46±| 


4 


J1223.3+7954 




6 


0.7/0.8/1.5 


61+e 


4 


J1347.0-2956 




6 


0.7/0.7/1.4 


55±| 


2 


J1410.4+7411 


18 


6 


0.4/0.9/1.3 


35 ±2 


4 


J2257.9-3646 




6 


0.2/0.7/0.9 


18+? 


4 


J2347.2+0707 


19 


5 


0.7/0.6/1.3 


84^° 


4 




J0143.6-5844 


30 


23 


0.5/1.3/1.8 


53 ±4 


4 


J0305.0-1602 


16 


6 


0.5/1.0/1.5 


39+? 


4 


J0338.2+1306 


38 


16 


2.2/1.0/3.2 




4 


J1410.4+7411 


33 


10 


0.8/1.6/2.4 


36 ±2 


2 



Table 5: Top: 24-month data: Number of 7-ray photons, listed for detected energy bins 
between 3-10 GeV and 10-300 GeV, within a radial region of 0.5° around the source's position. 
Additionally, the expected numbers of background photons from Galactic foreground (gal) 
and isotropic background (iso) between 10-300 GeV are given. The last columns list the 
energy of the photon with highest energy, and its corresponding event classification assigned 
by LAT data reconstruction (2: SOURCE, 4: ULTRACLEAN). Bottom: Same as above, 
for 42-month of data: Number of 7-ray photons in the 42-month data set for spectrally 
preselected source candidates, see section 3.3. 



for the spectrally selected list of candidates after 42 months. Additionally, the number of 
background photons N ^ expected within 0.5° around the source is given for the 10-300 GeV 
interval. These values were derived by fitting a Rol of 1° x 1° centered on the position of 
the Sol, fixing the normalizations (and power-law correction) of the Galactic foreground and 
isotropic background templates to the best-fit values obtained from the entire 0.1-300 GeV 
fit. 8 On average, approximately seven photons above 10 GeV have been detected after 24 
months, containing between one to two background photons within a radius of 0.5°. We note 
that only 2FGL J0338. 2+1306 has been detected above 100 GeV. 

The potential variability of the source flux was tested with an unbinned Kolmogorov- 
Smirnov (KS) test [53]. This test is already valid for a small number of photon counts, 
in distinction to the binned chi-square method used for the 2FGL catalog. The empirical 
cumulative distribution function of the arrival times of individual photons is compared to 
a uniform distribution, taking the (possibly varying) exposure into account. The resulting 

8 Owing to the quadratic-shaped Rol, the numbers were multiplied with 7r0.5 2 w 0.785 to correct for a 
circular Rol. 
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2FGL Name 






95% upper 


limit 0f fdeg] 


3-300 GeV 10- 


-300 GeV 


3-300 GeV 


10-300 GeV 


J0031.0+0724 




0.51 




0.53 


J0116. 6-6153 




0.17 




0.55 


J0143. 6-5844 


0.97 


0.83 


0.25 


0.38 


J0305. 0-1602 


0.006 


0.19 


0.45 


0.73 


J0312. 8+2013 




0.05 




0.50 


J0338. 2+1306 


0.15 


0.57 


0.48 


0.60 


J0438.0-7331 


0.28 


0.33 


0.30 


0.50 


J0737.5-8246 




1.00 




0.35 


J1223.3+7954 




0.23 




0.35 


J1347.0-2956 




0.40 




0.70 


J1410.4+7411 


0.78 


0.66 


0.38 


0.60 


J2257.9-3646 




0.08 




0.90 


J2347.2+0707 


0.62 


0.65 


0.28 


0.75 




J0143.6-5844 


0.89 


0.46 


0.17 


0.25 


J0305.0-1602 


0.004 


0.04 


0.45 


0.70 


J0338.2+1306 


0.02 


0.82 


0.23 


0.38 


J1410.4+7411 


0.47 


0.97 


0.33 


0.47 



Table 6: Probability for temporally constant 7-ray emission (left columns) and upper limit 
on the intrinsic angular extent of the signal (95% c.L, right columns). The quantities are 
listed for the 10-300 GeV band and, in case of a detection between 3-10 GeV, 3-300 GeV 
band. Quantities derived from the 24 (42)-month data are shown in the top (bottom) panel. 



probabilities -P C onst for the temporal photon distribution to be consistent with a constant flux 
are listed in table 6 for both the 3-300 and inclusive 10-300 GeV interval. In particular, we 
find 2FGL J0305. 0—1602 to show indications for variability (-P C onst = 4%o). 

We used a likelihood-ratio test to probe for intrinsic spatial extent. The corresponding 
likelihood function is given by L(8 S ) = —2J2iLi ^ n [Pdet( x i — x; 9 S ) + b], where pdet( x ; #s) fol- 
lows the probability distribution function (PDF) for a photon to be detected at x, x denotes 
the best-fit position of the Sol, and b denotes the (flat) PDF of the underlying background 

x — x for an intrinsically extended 7-ray emitter is the convolution 
of Ferrm-LAT's point spread function (PSF) ppsf (version P7_V6) with the intensity profile 
Pint of the emitter, p^et = Ppsf * Pint- The intensity profile p mt of a DM subhalo follows 
its line-of-sight integrated squared density profile. Similar to the approach followed in Pa- 
per I, the DM density profile of a subhalo was assumed to follow the spherically symmetric 
Navarro-Frenk- White (NFW) profile p(r) oc [r/r s (1 + r/r s ) 2 ] [54], where r denotes the 
distance to the center of the halo with a characteristic value at r s . Since 87.5% of the total 
luminosity is produced within r s , it serves as a convenient proxy for the intrinsic subhalo 
extent. For a subhalo at distance D, this corresponds to the angle 9 S ~ r s /D. We remark 
that 68% of the total luminosity are produced within #68 — 0.46 6 S , which is more convenient 
for comparison with observational data. 

The minimum L m ; n of the likelihood function L arises for the extension parameter 
6 S fitting the photon distribution best. Applying the theorem of Wilks [51], the quantity 
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AL = L — L m in follows a chi-square distribution with one degree of freedom, with additional 
terms of the order of 1/iV 1 / 2 , which are important for a small number of counts (see also 
[55]). 

We find no source candidate which shows indication for an intrinsic angular extent. 
Note that this result is in agreement with ref. [28], that also searched for angular extended 
sources. For each source, we present upper limits on 9 S at 95% confidence level in table 6. 
For the 24-month data between 10-300 GeV, note hereby that the given confidence level is 
not precisely defined, owing to the low number statistics which might affect the chi-square 
distribution. In general, the upper limits range from 0.2° to 0.9°, where the most constraining 
ones can obviously be derived from the largest data sets (42 months, 3-300 GeV). 

3.5 Multi-wavelength counterparts 

3.5.1 Catalog data 

Although the sources in table 1 are cataloged as unassociated, we carried out a dedicated 
counterpart search. 9 In particular, radio and X-ray sources positionally located inside the 
95 %-confidence-level uncertainty contour (listed in the 2FGL) of the 7-ray source have been 
searched for, providing counterpart candidates in case of a non-DM origin. The resulting 
counterpart candidates are listed in table 9 (see appendix B), ordered by increasing angular 
separation from the 2FGL position. For every radio or X-ray source, the USNO-B1.0 catalog 
[56] was searched for a matching optical counterpart. 

For all but one source (2FGL J1410. 4+7411), at least one radio source is located nearby. 
However, note that owing to the high sensitivity (therefore low confusion limit) of the NVSS 
[57] and SUMSS [58], which were used for radio associations, the high radio source density 
(~0.2arcmin -2 ) yields a large number of by-chance associations in the rather large positional 
uncertainty of the 7-ray sources (see table 1 and figure 2). 

3.5.2 WISE data 

The sources 2FGL J0312. 8+2013, 2FGL J0737.5-8246, and 2FGL J1347.0-2956 have re- 
cently been associated with blazar candidates selected from the WISE survey [59] on the 
basis of their mid-infrared (IR) colors [7]. At least one WISE object located in the positional 
uncertainty of these 7-ray sources has been found to fulfill the IR properties of (7-ray) blazars, 
i.e., lying in the WISE gamma-ray strip. We adopted this approach [7, 60, 61] to classify 
IR-counterpart candidates of 7-ray sources by their mid-infrared color-color properties. We 
focus on the IR counterparts of the four spectrally selected 2FGL candidates, together with 
the candidate selected in Paper I, 2FGL J0031. 0+0724. For each of these five 2FGL sources, 
an infrared counterpart candidate is listed in the WISE catalog [62], which is positionally 
coincident with the established radio and X-ray associations (see section 3.5.3, table 9, and 
section 3.6). The associations are listed in table 10 (appendix B), together with their infrared 
magnitudes W1,W2,W3, and WA. As described in [60, 61] and references therein, the pop- 
ulation of blazars spans a distinct region in the infrared color space of WISE objects, the 
WISE Blazar Strip. In particular, the sample of known 7-ray emitting blazars populates the 
WISE gamma-ray strip (WGS), which is a subspace of the WISE Blazar Strip. The WGS can 
be parametrized in two (overlapping) subregions, one which is dominantly populated by 7-ray 

9 Apart from preselected catalogs, the archives NASA/IPAC Extragalactic Database (NED, 
http:/ /ned. ipac.caltech.edu/) and HEASARC (http:/ /heasarc. gsfc.nasa.gov/cgi-bin/W3Browse/w3browse. pi) 
have been queried. 
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Figure 1: (a): Infrared color-color diagram of the WISE associations in table 10 (blue 
circles). The WGS subregion of BZBs is bordered with the dashed line, while the reference 
sources used in [61] are indicated with the gray diamonds, (b): Same as (a), comparing to 
the WGS subregion of BZQs. 



emitting BL Lacs (BZBs), while the other one is dominated by flat spectrum radio quasars 
(BZQs). In figure 1, we compare the infrared color-color diagram (W2 — W3,W1 — W2) of 
the WISE associations in table 10 with the BZB and BZQ regions of [61]. 10 We find that the 
WISE associations for all five 2FGL sources hint for a BL Lac origin of the 7-ray emission, 
as being consistent with the BZB region of the WGS. 

3.5.3 Swift-UVOT /X.RT data 

Besides the available catalogs of known X-ray sources, we have searched for unpublished 
archival observations of the remaining four objects. In a dedicated campaign [63], UV and 
X-ray follow-up observations of unidentified 2FGL sources have been carried out with the 
UVOT and X-ray telescope (XRT, 0.2-10keV) aboard the Swift satellite [64, 65]. 

The photometric UVOT data have been extracted from the products of the standard 
pipeline using the HEAsoft 6.12 software package in combination with the calibration files 
(2012-04-02). The in-orbit calibration procedures are described in detail in [66] . The standard 
aperture of 5 arcsec has been used for all filters to extract the background subtracted flux 
with uvotsource. 

For the corresponding XRT data sets (see table 7 for details), calibration and screening 
(xrtpipeline) of the data acquired in photon-counting (PC) mode was done using standard 
screening criteria, along with the current release of calibration files (2012-04-02). Data were 
reduced with the HEAsoft 6.11 software package. We used Ximage for source detection and 
Xspec (version 12.7.0) for spectral fitting. The probability limit for a background fluctuation 
was set to the 5<7-level and we required a signal-to-noise ratio of S/N > 4. Positions and 
corresponding uncertainties were derived with xrtcentroid. For each source position, ancillary 
response functions needed for spectral fitting were derived with xrtmkarf incorporating PSF 
correction. The circular on-source regions contained about 90% of the PSF (corresponding 



Since most of the WISE sources tabulated in 10 have not been detected in the WA band, we did not 
attempt to use the color-color projections including W4. 



FoV 



Obs. ID Obs. year Exposure [ks] 



2FGL J0143.6-5844 
2FGL J0305. 0-1602 
2FGL J0338.2+1306 
2FGL J1410.4+7411 



41274 2010 4.4 

41286 2011 3.2 

41292 2010 4.1 

47219 2012 3.5 



Table 7: Archival Swift-'KKT data for the celestial regions of the preselected source candi- 
dates. The columns list the observation ID, the observation year, and the total exposure in 
ks. 



to a radius of ~47"), while the background was derived from appropriate off-source regions 
with radii between 3' and 5'. All X-ray sources were spectrally fit with a power-law model 
corrected for photoelectric absorption. The hydrogen column density Ah was fixed to the 
nominal Galactic value, calculated from the Leiden/ Argentine/Bonn (LAB) HI survey [67]. 
For faint sources (S/N < 6), the power-law index T was fixed to 2.0. To achieve sufficient 
fit quality, the spectral channels were grouped, requiring a minimum of 5 (for S/N < 15) or 
10 counts per bin, respectively. We used the Cash-statistic for spectral fitting, to properly 
treat the low count statistic. 

In general, we find new X-ray sources in every Swift field of view (FoV), detected 
with comparably hard spectra having (unabsorbed) fluxes at the level of 10 -13 erg cm -2 s -1 
between 0.2 and 2keV. For all selected candidates but 2FGL J1410. 4+7411, some X-ray 
sources are positionally consistent with the cataloged 7-ray uncertainty, and additionally 
with the radio detections mentioned above (see table 9 in appendix B). Therefore, they 
provide convincing counterpart candidates for the 7-ray sources, cf. Paper I. For reference, 
positional and spectral parameters of every X-ray detection are listed in table 11 (appendix 



3.6 Discussion of preselected candidates 

Below, we provide a discussion of every preselected source candidate, based upon the results 
obtained in sections 3.3 to 3.5. Additionally, the updated results on 2FGL J0031. 0+0724 are 
summarized. 

Positional and spectral features of the four selected 2FGL candidates and their corre- 
sponding multi-wavelength associations (established in section 3.5) are summarized in figures 
2 and 3, respectively. In figure 2, the 2FGL as well as 42- month best-fit position and their 
corresponding uncertainties are overlayed over the photon image measured with Swift-~KRT , 
and plotted together with the positions of radio and X-ray sources. For all four 2FGL sources, 
the updated best-fit position shifts by a few arcmins, where the largest shift was found for 
2FGL J0338. 2+1306 (5.2'). Figure 3 compares the multi-wavelength data with the aver- 
age spectral energy distribution (SED) of a high-energy peaked blazar (HBL), which has 
been adapted from [72, 73] for particular redshifts z. Note in this context that the multi- 
wavelength data have not been taken contemporaneously. The redshift of each source was 
estimated from the distance modulus tor — Mr, where itir denotes the magnitude in the 
USNO-B1.0 R-band (table 9), and we assumed the detected optical emission to originate 
from a standard giant elliptical host galaxy with an absolute magnitude Mr = —23.1 [74]. 
We assumed a vanishing K-correction, i.e., a power-law spectrum with index a = T — 1 = 1. 
We emphasize that this method only provides a rough estimate under the given assumptions, 



B). 
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Figure 2: Best-fit position and uncertainty contour of (a) 2FGL J0143.6— 5844, (b) 
2FGL J0305. 0-1602, (c) 2FGL J0338. 2+1306, and (d) 2FGL J1410.4+7411 for 24 and 42 
months of Fermi-LAT data. The cataloged position is marked with the red "+", while the 
dashed red line borders its uncertainty ellipse (95% c.l.). The black "x" marks the 42-month 
position, and the solid black line its uncertainty contour (95% c.l.). An image of X-ray pho- 
tons (Swift-XRT), smoothed with a Gaussian (7"), is shown in the background. Note that 
the region of (c) has not been entirely observed with Swift-XRT, and different z-axis scales 
are used to improve readability (i.e., (a) log, (b) sqrt, (c) linear, and (d) linear). Positions 
of radio sources (NVSS) are indicated with dark-green boxes, discovered X-ray sources with 
blue arrows. Note that the boxes's size does not reflect the positional uncertainty. 



while a precise determination of z requires spectroscopic data in the optical band. In the 
very high-energy regime, emitted photons are absorbed through 7-7 pair production, which 
was calculated using the EBL model provided in [75] . 

2FGL J0031. 0+0724. The 24-month Fermi-LAT data of this source has been intensively 
studied in Paper I. The analysis carried out here does not demonstrate a preference for an 
exponential cutoff in the 42-month data set. With a significance of ~ 2a, the lightcurve 
of high-energy photons (10-300 GeV) is consistent with a temporally variable source, and 
no indication for angular extended emission was found. In Paper I, we already claimed a 
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Figure 3: Spectral energy distribution (SED) of (a) 2FGL J0143. 6-5844, (b) 
2FGL J0305.0-1602, (c) 2FGL J0338.2+1306, and (d) 2FGL J1410.4+7411, assuming 
the multi-wavelength associations discussed in the text. Included multi-wavelength data, 
from low to high frequency: radio (NVSS, 1.4 GHz; black triangle), infrared (WISE, 
Wi,W3,W2,Wl; red diamonds), optical (USNO-B1.0, R, B; green triangles), ultra-violet 
(GALEX, NUV, FUV; violet boxes; see http://galex.stsci.edu/GR6/; Swift-\JVOT, U, 
UVW1, UVM2, UVW2; dark-golden points), X-ray (Swift, 0.2-2keV; blue line), 7-ray 
(Fermi-LAT 2FGL, 0.1 — lOOGeV; red line and circles). The optical and UV data have been 
dereddened using E(B — V) from [68] and assuming Ry = 3.1 (see [69] for details). Arrows 
indicate upper limits (95% c.l.). Statistical uncertainties of the X-ray and 7-ray spectra are 
indicated by the corresponding shaded areas [70]. The orange line shows the sensitivity of 
the planned CTA observatory for 50 hours of observation [71]. For comparison, the solid 
black line shows the average SED of a high-frequency peaked blazar (HBL), adapted for the 
estimated redshifts z. The HBL-SED is normalized to the radio flux, and the energy flux vf v 
is plotted in the frame of a potential observer. The HBL-SED has been corrected for EBL 
absorption (see text for details), while the dotted black line shows the SED for a vanishing 
EBL. 
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possible association of the source with a faint radio source (12mJy), positionally coincident 
with a faint X-ray source, / unabs (0.2— 2 keV) ~ 2.1 x 1CP 13 erg cm -2 s _1 . The spectral energy 
distribution (SED) is consistent with a BL Lac origin (see Paper I), which is supported by 
the photometric infrared data of WISE J003119. 70+072453. 6. 

2FGL J0143.6— 5844. The 42-month data indicates the 7-ray spectrum to be preferen- 
tially fit with an exponential cutoff with ~ 3a. Furthermore, the lightcurve is preferred to 
be steady, but the 7-ray emission is consistent with a point-source. Multi-wavelength obser- 
vations show the source to be associated with a 27mJy radio source (SUMSS), positionally 
coincident with the bright X-ray source SWIFT J014347.3-584551, / unabs (0.2 - 2 keV) ss 
1.4 x 10~ n erg cm~ 2 s . Together with the simultaneously measured X-ray flux, the UV 
data seems to indicate variable emission when compared with the optical data. The infrared 
data supports a BL Lac scenario, which is also consistent with the multi- wavelength SED, 
see figure 3a. 

2FGL J0305. 0—1602. While the 7-ray data shows initial indication (at the 2.5a level) 
for a spectral cutoff at ~ 10 GeV, the temporal photon distribution of the source excludes 
a constant flux with ~ 99% confidence. Multi-wavelength searches indicate the source to be 
positionally associated with PKS J0305— 1608, showing a radio flux at the Jy level (NVSS) 
and an X-ray flux of / unabs (0.2 — 2 keV) rj 2.3 x 10~ 13 ergcm~ 2 s _1 . The photometric in- 
frared data is sparsely consistent with a BL Lac scenario, while the high radio flux is not in 
accordance with the expectation from a high-energy peaked BL Lac (see figure 3b). 

2FGL J0338. 2+1306. The source was initially selected based upon the 24-month data set, 
preferring a spectral cutoff with a significance of ~3c As shown in figure 2c, gaining photon 
statistics revealed a large positional shift fitting the 42-month data, and the initial indication 
for a spectral cutoff vanished. The updated data set also indicates a variable 7-ray flux at 
the 2a level. The improved positional accuracy allows to associate 2FGL J0338. 2+1306 with 
a radio source (15.1 mJy), which is positionally coincident with SWIFT J033829. 2+130217, 
see figure 2c. Its WISE counterpart suggests a BL Lac origin, in accordance with the entire 
multi- wavelength emission (see figure 3c). 

2FGL J1410.4+7411. The 42-month 7-ray data of this source prefers its spectrum to 
be fit with a power-law (r = 1.5) with exponential cutoff, with a significance of ~ 3. 3d. 
The lightcurve is consistent with steady emission, but no indication for an angular extent 
was detected. While also being counterpartless after 24 months, the position computed 
from the larger photon sample seems to shift towards the very faint X-ray source SWIFT 
J141234. 8+741153, see figure 2d, / unabs (0.2-2 keV) « 6.3 x 10~ 14 erg cm" 2 s" 1 . 11 Note that 
the Swift source has no radio counterpart, which might reflect the lacking sensitivity of radio 
surveys. Assuming the Swift source to be the correct X-ray association, its infrared counter- 
part WISE J141235. 75+741158.0 would indicate a BL Lac origin. Additionally, compared to 
the USNO and GALEX data, the Swift-WOT observations indicate variable emission. 

4 Discussion and conclusions 

In this work, we investigated the unassociated 7-ray source population of the Fermi-LAT 
second year point-source catalog for sources potentially originating from DM subhalos. Basic 

"Note that the X-ray source XMMSL1 J141002.6+740744 (/ unabs (0.2-2 keV) w 2 x 1(T 12 erg cm" 2 s" 1 ) [76] 
is located just outside the 2FGL 95% uncertainty contour, but does not appear in the SWIFT observations. 
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catalog selection revealed 13 high-latitude sources, with hard 7-ray spectra detected above 
10 GeV, and lacking indication for temporal variability. Using 3.5 years (42 months) of Fermi- 
LAT data, we developed a statistical test to probe the candidates for spectral consistency with 
self-annihilating DM (a power-law with exponential cutoff) . The high-energy spectra of a sub- 
set of 4 sources were found to be preferentially fit by power-laws with exponential cutoff (i.e., 
2FGL J0143.6-5844, 2FGL J0305. 0-1602, 2FGL J0338. 2+1306, and 2FGL J1410.4+7411), 
with significances between 2.5a and 3.3o\ All sources were tested for temporally constant 
and spatially extended 7-ray emission. The 7-ray emission of 2FGL J0305. 0—1602 shows 
a ~ 3<7 indication to be temporally variable, while none of the 13 sources shows indica- 
tions for angular extended emission. Multi-wavelength studies were conducted to search 
for associations, using refined positional information based upon 3.5 years of Fermi-LAT 
data. For 2FGL J0143. 6-5844, 2FGL J0305.0-1602, and 2FGL J0338. 2+1306, we estab- 
lished clear associations detected in the radio, infrared, optical, UV, and X-ray bands, while 
2FGL J1410. 4+7411 is indicated to be associated to a faint X-ray source, which has also 
been detected in the infrared, optical, and UV band. 

With the exception of 2FGL J0305. 0—1602, the infrared color-color data of the three 
other sources is fully consistent with the population of BL Lacs detected with Fermi-LAT . 
In addition, such a scenario would be favored by the multi-wavelength data, in particular 
by the faint radio, hard X-ray, and hard 7-ray emission, indicating a scenario of a high- 
frequency peaked BL Lac (see figure 3 and Paper I for details). For all three cases, we note 
that the 7-ray flux predicted in the average model is below the Fermi-LAT measurement. 
This might indicate a sample bias, meaning that the sources have been detected in a flaring 
state. Within the errors, the infrared association of 2FGL J0305. 0—1602 might also indicate 
a BL Lac origin, while in particular its bright radio counterpart may point towards a differ- 
ent scenario. Finally, the recent study in [8] has attempted to classify the entire sample of 
unassociated Fermi-LAT sources, distinguishing between AGN-like and pulsar-like sources. 
Using a Random Forest (RF) classifier trained on cataloged 7-ray properties of associated 
Fermi-LAT sources, the major fraction of the unassociated sources are predicted to be AGN, 
and no significant outliers have been found. All 13 candidates we selected in table 1 are 
suggested to be AGN. Consistently, the recent RF classification in [9] assigns all of them to 
originate from BL Lacs. In particular, this strongly supports the BL Lac origin of the four 
spectrally selected candidates. Note, however, that since we predict 7-ray sources originating 
from DM subhalos to be particularly faint, their cataloged spectral and localization param- 
eters suffer from large statistical uncertainties. Therefore, DM subhalos might hide in the 
sample of sources classified by the RF algorithm, emphasizing the necessity of the in-depth 
investigations carried out in this paper. 

In conclusion, we find no unassociated 7-ray source in the 2FGL catalog which is favored 
to originate from a subhalo driven by self-annihilating WIMPs at a mass scale above 100 GeV. 
However, we conclude that, among all candidates, the source 2FGL J1410. 4+7411 would 
be the most interesting, owing to 7-ray properties which might prefer a DM origin and 
a high uncertainty about its association. From the final source sample, we can exclude 
2FGL J0305. 0—1602 (owing to variable 7-ray emission). We note that no source can be 
firmly excluded by spectral properties. We find the remaining candidates to most likely 
originate from faint BL Lacs. 

Albeit our prediction of a BL Lac origin, we investigated all source candidates in the 
context of the recently claimed evidence for a line-like feature at ~ 130 GeV in the Galactic 
Center region [38, 77-83]. In case of a self-annihilating DM scenario, DM subhalos will also 
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appear with a 7-ray line at ~ 130 GeV, and searches have been started in [80, 84-86]. Except 
for 2FGL J0338. 2+1306, having a photon at 152 GeV, we note that none of the other sources 
has been detected above 100 GeV. 

The study presented here has clearly outlined the problems of identifying faint Fermi 
sources. Difficulties are mainly related to the limiting collection area of the LAT at the high 
energy end, resulting in a small number of photons. This implies the consequently large 
uncertainty of source locations, which in turn leads to source confusion at the faint end of 
source populations. Likewise, the small number of photons limits the ability to determine 
spectral and temporal properties at the high energy end of the Fermi-LAT response with 
sufficient accuracy to distinguish source models. Finally, from the theoretical point of view 
in the considered DM scenario, the entirely limited number of detectable subhalos prohibits 
conclusive population studies. 

Focussing on high-latitude sources with fluxes at the level of the studied ones, we em- 
phasize that at least some of these issues can be addressed with larger data sets based upon 
longer observations. The correspondingly large signal-to-background ratios allow us to im- 
prove the positional accuracy by a factor of ~2 with 10 years of Fermi-LAT data, for instance, 
and to monitor the temporal photon distribution over longer time periods. The improvement 
in sensitivity might push the number of detectable subhalos to 0(5). 

In particular, the issues can also be addressed with pointed follow-up observations in the 
very high-energy (VHE) band. The large collection areas provided by imaging air Cherenkov 
telescopes (IACTs) for energies above ~ 50 GeV, such as H.E.S.S.-II [87, 88], MAGIC stereo 
[89, 90], VERITAS [91, 92], and, in particular, the planned Cherenkov telescope array (CTA) 
[71, 93, 94], allows the detection of a larger photon sample, significantly reducing positional 
and spectral uncertainties. If detected in the VHE, source candidates therefore can be easily 
associated or even identified, see [27]. 

As a final remark, a potential successor of Fermi-LAT such as GAMMA-400 [95, 96] 
will significantly improve the observable energy range (100 MeV — 3 TeV), angular resolution 
(~0.01° at 100 GeV), and energy resolution (~1% at 100 GeV). The launch of GAMMA-400 
is planned for 2018. For the case of unidentified Fermi-LAT sources, such a telescope will 
constrain their celestial position with enhanced precision. 
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A Probability distribution of TS exp 

To determine the probability density distribution (pdf) of TS cxp (eq. 3.1) in the null- 
hypothesis Ho, we used bootstrap Monte Carlo simulations [53]. Assuming a pure power-law 
spectrum (i.e., Hq), we simulated data of the Rol corresponding to 2FGL J0338. 2+1306, 
which was exemplarily selected to check whether an exponential cutoff is preferred by the 
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actual data set (see table 3). The 25 000 simulated data sets were then analyzed with the 
framework described in section 3.2 and 3.3, therefore calculating TS exp in exactly the same 
way as in the actual data analysis. Each individual analysis procedure was checked for non- 
converging behavior (~14%), and the final data set of simulations was cleaned accordingly. 

The simulations were done in a two-step approach: First, we used gtobssim to simulate 
five individual 24-month data sets of the Rol centered on 2FGL J0338. 2+1306 (between 
100 MeV and 300 GeV). All sources in the Rol were modeled with the cataloged power-law 
spectra, since the current version of gtobssim does not accept log-parabola spectra (which 
are sometimes preferred). Background models were implemented as described in section 3.2, 
and we used the actual spacecraft file of the first 24 months of data-taking. We checked the 
consistency between the analysis results from the simulated data and the actual data set. 
As a second step, the five individual simulations were merged to one data set. Applying the 
bootstrap technique to the merged data, we generated 25 000 individual 24-month data sets 
to be used in our analysis. Again, we checked that the analysis reproduced the actual data 
well. 

Figure 4a shows the pdf of TS cxp , fixing the index T of the alternative hypothesis Hi 
(power-law with exponential cutoff) to T = 1.5 and T = 0.35, respectively. The simulation 
shows that the pdf does indeed not follow a ^/^-distribution [52] for both negative as well 
positive TS ex p values. Rather, we find asymmetric pdfs with maxima and large tails at 
negative TS exp values. 12 These features are more pronounced for T = 0.35, owing to the 
sharply peaked 7-ray spectrum. For the positive half-plane, the p-values p{> TS eX p) (the 
significance of the test statistic) are shown in figure 4b. The p-values are compared to 
the prediction of Chernoff's theorem [52], where the test statistic in the positive half-plane 
should follow 0.5 [<5(TS exp ) + Xi (TSexp)] , with 5(x) the delta-funct ion and Xii x ) the chi- 
square distribution with one degree of freedom. For T = 0.35, we find that for large, positive 
TS exp the pdf approximately follows Chernoff's theorem, while the pdf of a T = 1.5 spectrum 
does not. The resulting significances corresponding to selected TS exp values are shown in 
table 8. For the index T = 1.5 (0.35), a significance of 2a corresponds to TS cxp = —6 (—20) 
and TSexp = 2(2), respectively, while the 3a contour is given by TS 0X p = —25 (—35) and 
TSexp = 4(7). 



12 The maxima are at TS cxp = 



-0.25 (-2.25) for F = 1.5 (0.35). 




Figure 4: (a): Probability density distribution of the test statistic TS ex p in the null hypoth- 
esis Hq. The power-law index of the alternative hypothesis H\ was fixed to T = 1.5 (top) and 
T = 0.35 (bottom). The dotted vertical line indicates TS exp = 0. (b): P-value p(> TS exp ) for 
the positive half-plane of the distribution, starting from its corresponding maximum. The 
curves for T = 1.5 (blue line) and F = 0.35 (red line) are compared to a Xi/2 distribution 
(green line). 
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Table 8: Probability (p-value) to randomly find the test statistic to be less (left table) or 
larger (right table) than a certain TS exp value, given in Gaussian sigma. The p-value is listed 
for both assuming the index T of the power law with exponential cutoff to be 1.5 and 0.35, 
respectively. 
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